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Abstract. Considering a plasma with an initially weak large scale field subject to nonhelical turbulent stirring, Zeldovich 

(1957), for two-dimensions, followed by others for three dimensions, have presented formulae of the form (b 2 ) = }{Rm)B 2 . 
Such "Zeldovich relations" have sometimes been interpreted to provide steady-state relations between the energy associated 
with the fluctuating magnetic field and that associated with a large scale or mean field multiplied by a function / that depends 
on spatial dimension and a magnetic Reynolds number Rm- Here we dissect the origin of these relations and pinpoint pitfalls 
that show why they are inapplicable to realistic, dynamical MHD turbulence and that they disagree with many numerical 
simulations. For 2-D, we show that when the total magnetic field is determined by a vector potential, the standard Zeldovich 
relation applies only transiently, characterizing a maximum possible value that the field energy can reach before necessarily 
decaying, in relation to a seed value B. In 3-D, we show that the standard Zeldovich relations are derived by balancing 
subdominant terms. In contrast, balancing the dominant terms shows that the fluctuating field can grow to a value independent 
of Rm and the initially imposed B, as seen in numerical simulations. We also emphasize that these Zeldovich relations of 
nonhelical turbulence imply nothing about the amount mean field growth in a helical dynamo. In short, by re-analyzing 
the origin of the Zeldovich relations, we highlight that they are inapplicable to realistic steady-states of large Rm MHD 
turbulence. 
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1. Introduction 

Stars, galactic interstellar media, the hot plasma of galaxy 
clusters, and accretion disks all contain magnetohydrody- 
namic (MHD) turbulence. Over the past 50 years, the ubiq- 
uity of MHD turbulence and magnetic fields in astrophysics 
has stimulated an effort to understand the resulting magnetic 
and kinetic energy spectra (e.g. Biskamp 1997; 2003; Bran- 
denburg & Subramanian 2005). 

Simplified semi-analytic models combined with idealized 
numerical experiments with restricted boundary conditions 
have been used in combination, to help bridge the gap be- 
tween ignorance and realistic astrophysical MHD turbulent 
settings. One such class of studies focuses on the 3-D ampli- 
fication of a weak seed magnetic field by a nonhelical turbu- 
lent flow in a closed volume or periodic box (e.g. Kazanstev 
1968, Kulsrud & Anderson 1992; Maron & Blackman 2002; 
Haugen, Brandenburg, Dobler 2003; Haugen, Brandenburg, 
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Dobler 2004; Maron, Cowley, Mc Williams 2004; Schekochi- 
hin et al. 2004). Starting with an initial seed field, the sys- 
tem is forced with an incompressible turbulent flow. The 
key questions are: how much magnetic energy grows and 
what determines its saturated value and spectrum? Two initial 
value problems can be distinguished: (1) An initially weak 
(magnetic energy <C turbulent forcing energy) seed magnetic 
field with a finite power at wavenumber k — 0, subject to 
forced, non-helical isotropic turbulence at forcing wavenum- 
ber ki ~ 5fci (e.g. Haugen & Brandenburg 2004) (2) An ini- 
tially weak seed random magnetic field with no k = field, 
subject to forced, non-helical isotropic turbulence at fcf ~ 5 
(e.g. Maron et al. 2004; Schekochihin et al. 2004). 

The two initial value problems above address the 3-D 
small scale dynamo, defined as field amplification at or below 
the forcing scale. Simulations show that for nonhelical ran- 
domly forced flows, no significant field amplification above 
the scale of the forcing occurs (e.g. Maron & Blackman 2002; 
Haugen et al. 2003, 2004; Maron et al. 2004). When a k = 
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field is imposed and a numerical simulation is performed with 
periodic box, this mean field cannot change due to the bound- 
ary conditions. If the turbulence is instead forced helically at 
k = fcf, the k = field still cannot grow, but the field up 
to k = 1 can and does grow (e.g. Pouquet, Leorat & Frisch 
1976; Brandenburg 2001; Maron & Blackman 2002; Black- 
man & Field 2002). (The explicit difference between the he- 
lical and non-helical forcing can be seen in the set of simu- 
lations by Maron & Blackman (2002).) Any significant field 
amplification on scales above the forcing scale, involves the 
large scale dynamo mechanism, and requires helical forcing. 
The nonlinear helical dynamo is the dynamical generaliza- 
tion of the original large scale dynamo (e.g. Moffatt 1978; 
Parker 1979) proposed to model the large scale fields of ro- 
tating stratified systems like the sun or Galaxy. 

The present focus is on the nonhelical dynamo saturation. 
Two key results emerging from 3-D numerical simulations of 
the nonhelical dynamo are that ( 1 ) no significant field ampli- 
fication occurs at k < k{ and (2) the total turbulent magnetic 
energy integrated for k > kf grows to within a factor of a few 
of equipartition with the incompressible turbulent kinetic en- 
ergy (when the magnetic Prandtl number Ptm =v/\ (where 
v is the viscosity and A is the magnetic diffusivity) is not too 
small; Haugen et al. 2004; Boldyrev & Cattaneo 2004. This 
is independent of whether there is an initial k = 0, or other 
large scale (0 < k < fcf) field present, as long as that initial 
large scale field is weak compared to the turbulent forcing 
(e.g. Haugen and Brandenburg 2004). In fact, the saturated 
end state of the nonhelical MHD turbulent dynamo spectrum 
is essentially independent of the shape of the seed spectrum 
(Maron et al. 2004). 

The approximate equipartition and its independence of 
an initially weak large scale field seem physically reason- 
able, but contradict a classic relation in MHD referred to as 
the "Zeldovich relation." This refers to the Zeldovich (1957, 
hereafter Z57) result derived for 2-D, and subsequently gen- 
eralized for 3-D using different approaches (e.g. Steenbeck 
& Krause 1969; Low 1972; Parker 1973; Krause & Radler 
1980; Zeldovich, Ruzmaikin, Sokoloff 1983 (hereafter Z8 3)). 
The Zeldovich relation is given by 

<b 2 } = /(i? M )B 2 , (1) 

where the magnetic field is B = b + B, with b the fluctuating 
field and B the mean or large scale field. Eq. Q is sometimes 
used to purport a steady state relation for the fluctuating (~ 
total) magnetic energy and the magnetic energy in B when 
the latter provides an initial seed. (We will see later that the 
distinction between B as a large scale field derivable from a 
vector potential vs. B as a mean fc = field has a particular 
implication for 2-D.) 

The function J{Rm) of the magnetic Reynolds number 
Rm is oc RTj for two dimensions (Z57; Moffatt 1978). In 
three dimensions, different values have been proposed ana- 
lytically: As discussed in Krause and Radler (1980, hereafter 
KR), for the large R M limit, Steenbeck & Krause (1969), 
Low (1972), Parker (1973) all find f(R M ) ~ R M , while 
/(Rm) ~ R\i for the low Rm limit. In contrast, in Z83, a 
more spectrally sensitive result is proposed. There /(-Rm) oc 



Iu(Rm) for three dimensions and Kolmogorov turbulence, 

5-3p 

and oc R M P m 3-D when the kinetic spectral index p ^ 5/3. 

Two issues immediately arise: First, why do the 3-D an- 
alytic calculations differ so strongly in the form of 
Second, having just discussed that modern simulations show 
that MHD turbulence saturates to a state in which the fluctu- 
ating field energy has a value independent of any large scale 
field, how can a relation like Q apply regardless of the spe- 
cific form of /(-Rm)? Where do relations of the form Q 
come from and how can we understand why they do not agree 
with simulated MHD steady states? We address these ques- 
tions in the present paper. 

Note that before the non-linear regime could be studied 
with numerical simulations, most dynamo theorists studied 
the kinematic regime. This refers to the regime in which the 
magnetic field amplified by a turbulent flow does not signifi- 
cantly back-react on this flow, so that the velocity is not sub- 
ject to strong Lorentz forces. This regime lends itself to an- 
alytic studies, but the limitations of the implications must be 
properly understood. The Zeldovich relations were derived 
from this perspective, but the specific reasons why different 
derivations appear to disagree and why none are applicable to 
the steady-state of fully developed MHD turbulence have not 
been succinctly elucidated. Doing so is our present goal. 

In section 2, we discuss the Zeldovich relation in 2-D, and 
explain that it cannot represent a steady state if the total field 
is derived from a vector potential. We show instead that it 
represents an absolute maximum that the fluctuating field can 
attain before necessarily decaying. In section 3, we show that 
in 3-D, the Zeldovich relations have been derived using ar- 
guments ignore key nonlinear growth terms and thus balance 
subdominant terms. We also show that, although the differ- 
ences in the 3-D forms of ([0 mentioned above can be at- 
tributed specifically do different approaches and inclusion of 
different subdominant terms, all of the analytical approaches 
drop the same dominant term. As a result, the generalized 
Zeldovich relation for 3-D MHD turbulence is inapplicable 
to steady-state large Rm MHD turbulence. We conclude in 
section 4. 

We emphasize throughout, that the Zeldovich relations 
were always derived in the context of nonhelical turbulence. 
They were therefore never meant to apply to a system in 
which B grows, and thus do not present any constraints on 
helical dynamo theory. This can lead to confusion because the 
Zeldovich relation sometimes appears in discussions whose 
main focus is on mean field helical dynamo amplification 
(e.g. Vainshtein & Cattaneo 1992). Our detailed discussion is 
therefore focused specifically on understanding the inappli- 
cability of the Zeldovich relation to steady-state nonhelical 
MHD turbulence. 

2. Understanding the 2-D Zeldovich Relation 

We first address the 2-D Zeldovich relation derived by Z57. 
Since no dynamo action can be sustained in 2-D, and the 
magnetic energy must ultimately decay, one immediately 
wonders what a relation between B and (b 2 ) could mean? 
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In fact Z57 is primarily a 2-D anti-dynamo theorem paper: 
we will see that the relation between B and (b 2 ) represents 
the maximum that (b 2 ) could possibly attain before it neces- 
sarily decays. 

We start with a basic argument showing that in 2-D, when 
the total magnetic field (mean + fluctuating) is determined by 
a vector potential, (i.e. B=VxA and no k — compo- 
nent for a periodic box) the total magnetic energy must decay 
when surface integrals vanish. We write the total the incom- 
pressible induction equation for the total magnetic field as 
(e.g. Moffatt 1978) 

d t B = Vx(vxB) + AV 2 B = B- Vv- v- VB + AV 2 B, (2) 

where v is the velocity flow, assumed to be imposed by exter- 
nal forcing with no mean component, and A is the magnetic 
diffusivity. 

Now following Z57 and restricting to 2-D (x, y) incom- 
pressible flow, we can then separate out the equation for the 
z-component 

aA = -v.V£ z +AV 2 B, (3) 

Multiplying by B z and spatially averaging for incompressible 
flows, we have 

d t (B 2 } = -d^Bl) + XBjidjB*) 2\((d k B z ) 2 ) 

= -2X((d k B z ) 2 ), (V 

where the latter equality follows when we assume surface in- 
tegrals vanish. Eq. @ shows that the b z contribution to the 
magnetic energy will decay. 

Because of the first term on the right of (|2ji, the x, y com- 
ponents of the magnetic field need not immediately decay, 
although they will eventually decay. To see this, note that the 
vector potential (defined such that B = V x A) satisfies 



d t A = v x (V x A) + AV 2 A - AV(V • A) - V<£. 



(5) 



As above, we can separate out the z-component, which deter- 
mines B x and B y , to obtain 

d t A z = -v-VA z + XV 2 A z , (6) 

since the z— derivatives vanish. Eq. has the same form 
as (|3}- Multiplying by A z and averaging, we again have for 
incompressible flow 

d t (A 2 ) = -di(viA 2 z ) + Xdi{diAl) - 2X((VA Z ) 2 ) = 



-2\{{VA Z ) 2 ), 



where the last equality again follows from ignoring the di- 
vergence terms for suitable boundary conditions. This shows 
that A z , like B z also ultimately decays. The ultimate decay 
of both B z and A z implies that the total magnetic energy 
also eventually decays. The asymptotic steady-state is one of 
negligible total magnetic energy, (B 2 ) ~ 0. This is the anti- 
dynamo theorem for 2-D (see also Moffatt 1978). 

It is important in this derivation that we took B = V x A 
because the first term on the right hand side of Q would 
not have emerged as a straightforward total divergence if we 
were considering a periodic box that included a strictly uni- 
form k = periodic field that has no well defined vector 
potential. (This latter case, was first studied numerically by 
Moss (1970)). 



Now let us take advantage of the positive nature of the 
last term of Eq.@ to rewrite that equation 



d t {A 2 ) = ~2X(A Z 2 )/S 2 (t), 



(8) 



where S(t) represents the dominant time dependent charac- 
teristic variation scale of A z . Allowing for the variation of 
this scale is important because although the magnetic energy 
ultimately decays in 2-D, random walk field line stretching 
can temporarily amplify the field. How much amplification 
takes place before the field decays? 

To answer this, consider that in 2-D, the field initially 
grows exponentially via the first term on the right hand side 
of J2), much like in 3-D (Kazanstev 1968; Parker 1979, Kul- 
srud & Anderson 1992). Early in the kinematic regime, 

B(t) ~B exp(</r), (9) 

where Bq is the root mean square of the initial magnetic field 
in the x—y plane, and r, for sufficiently steep kinetic energy 
spectra, is the correlation time of the forcing scale. It should 
be noted that 2-D forced turbulence exhibits a steeper kinetic 
energy spectrum than 3-D turbulence because of the tendency 
for enstrophy to inverse cascade (e.g. Davidson 2004). Ap- 
proximate inertial range energy spectra E(k) oc k~ 3 pro- 
posed by Batchelor (1969) Kraichnan (1970) seem to be 
consistent with 2-D simulations (Gotoh 1998; Lindborg & 
Alvelius 2000) though perhaps for slightly different reasons 
than originally thought (Davidson 2004). In 3-D, the kinetic 
energy spectrum E(k) oc fc~ 5 / 3 is such that the the r appear- 
ing in J32i would instead be closer to that of the smallest eddy 
scale than the forcing scale (Kulsrud & Anderson 1992)). 

As the field lines stretch from the turbulent motions in 
a confined area, an increase in field energy means that the 
field must build up on smaller and smaller scales, so S(t) de- 
creases. For incompressible flow, the field strength increases 
linearly with the length of the field line. The length of the 
field line varies inversely with the scale of the field variation. 
We can therefore write 

S(t) = L(B /B(t)) q = LExp[-qt/T], (10) 

where q is an index to allow deviations from a linear rela- 
tion, L the scale of the initial field Bq, and we have used (|9jl- 
Combining (jSJi and dlOt gives 



d t ln(A 2 z ) = -2XExp[2qt/r]/L 2 



(11) 



(7) so 



(12) 



ln((A 2 z )/(A 2 z ) ) = -(rX/qL 2 )(Exp[2qt/r} - 1) 
where we have defined the magnetic Reynolds number 

T 2 

Rm,l = xp> an d where the similarity follows for 2qt > r. 



Thus 

(A 2 Z ) = (Al) Q Exp[-{l/qR MtL )(B{t)/B Q ) 2 *]. (13) 

The field grows by stretching until the argument of Jl 31 be- 
comes > 1. This occurs when 

B{t) ~ B (qR M ,L) 1/2q - 

For the case in which Bq 
(6 2 ) 1/2 so (H3) then gives 



(b 2 ) 



2 - B\qR M ,L) 1/q , 



(14) 

B, we can also write B(t) ~ 
(15) 
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3. 1 Deriving the 3-D Zeldovich relation 



which for q — 1 is the 2-D Zeldovich relation. 
Using (1101 , Eq. (I14> also implies 
L 



(16) 



m ~ {qRKLL) 1 ' 2 ' 
This is the minimum scale at which the field energy could 
peak in the kinematic regime and at which point the vector 
potential and the total magnetic energy would rapidly decay. 
If the peak is away from this minimum scale, then the decay 
rate of A z is reduced: Using d!6i we can rewrite Jl 3I > as 



{A\) = (A 2 z ) Exp[-(S m /S(t)) 2 



(17) 



At the time Eq. Jl 5I > is satisfied the vector potential rapidly 
decays from Jl 3I > or J17I . and so does the total field energy. 
The Zeldovich relation J15l > therefore, simply provides an es- 
timate of the maximum value that the fluctuating field energy 
could obtain if it were to grow kinematically to this maximum 
(see also Pecko ver & Weiss (1 978)1. There is no steady state 
for which dl5t applies. Moreover, J15I is a kinematic result 
relevant for (v 2 ) > Rm,lB$. If instead (v 2 ) < Rm.lB 2 , the 
dynamical amplification of (b 2 ) is ultimately limited by near 
equipartition with the kinetic energy density well before the 
value implied by H5\ is reached. This latter statement also 
applies even when B represents a k = component; while 
the above proof of the total magnetic energy decay would 
not be valid in this case, equipartition with the kinetic energy 
would determine the limiting magnetic energy. 

Note also that the mean field enters ill 5b only if the initial 
field energy is determined by B. The same calculation would 
go through even if there were no initial mean field, in which 
case the initial field Bq would not be related to B. 

3. Understanding the 3-D Zeldovich Relation 

In 3-D, unlike 2-D, sustained dynamo action is allowed when 
the system is forced with a velocity flow, whether or not there 
is a k = mean field. As discussed in the introduction, sim- 
ulations of incompressible MHD turbulence show that the 
magnetic field typically grows to approach and saturate near 
equipartition with the turbulent kinetic energy in the steady 
state, independent of any weak initial large scale or k — 
mean field. In this section we show why the generalized Zel- 
dovich relation for 3-D does not account for this steady-state 
correctly. The distinction between large scale < k < kf 
vs. mean k = is not as essential for 3-D as it was in 2-D in 
what follows. 

3.1. Deriving the 3-D Zeldovich relation 

3.1.1. The steady-state Z83 approach 

We first follow Z83 and employ a spectral formalism where 
the kinetic energy spectrum E(k) is defined by 



E{k)dk = {v 2 )/2 
and satisfies 

2E(k) = v{kf/k = 



v 2 (h) (k 



(18) 



(19) 



where p is the spectral index, fc; is the wavenumber of the tur- 
bulent forcing scale, and v(k) is the turbulent speed at wave 
number k. Assuming a constant energy transfer rate gives 

v(k) 2 /r(k) = constant, (20) 

so the energy transfer time r(k) satisfies 

T (k)=r(h)(j) . (21) 
The magnetic energy spectrum M(k), is defined by 



M(k)dk = E M 



(22) 



where k\ is the resistive wave number and Em is the mag- 
netic energy density. 

Using B = B+b, where B is a fixed large scale field, and 
subtracting the mean of the magnetic induction equation (|2j 
from Eq. (|2|l gives the following equation for the fluctuating 
component of the magnetic field: 

d t b = V x (v x B) + V x (v x b) - V x (v x b) + AV 2 b.(23) 
We dot this equation with b and average to find 
1 



5 W 



1 'b-Vx(vxB)) + {b-Vx(vxb))+A(b-V 2 b).(24) 

Now we assume that the left hand side of ( 1241 is zero, thus 
assuming that a steady-state is maintained by the terms on 
the right side. Let us ignore the last term and non-local in- 
teractions between v and b, and follow the Z83 estimates 
of approximate magnitudes for the 1st and second terms on 
the right. To order of magnitude, the narrow band integrated 
Fourier spectrum for the first term on the right is given by 

(b- Vx (vxB)) (fe) ~b{k)kv(k)B. (25) 

Z83 replaces the narrow band integrated spectrum of the third 
term of J24b with a turbulent diffusivity, that is 



(b.Vx(vxb)) (fe) = - t / T (fc)fc 2 6 2 (fc), 
where 

v T {k) 



(26) 



(27) 



r(k)k 2 ' 

Assuming a steady state, setting Eq. ( I25> equal to minus 
Eq. ( I26> gives 

v(k)B = u T (k)kb(k). (28) 
Squaring this gives 



—2 



b 2 {k) = 2kM(k) 



v 2 {k)B _ 2E(k)B 
v^(k)k 2 v\ 



v\(k)k 



(29) 



Then using J 1 91 . J2 It and J27i we have 



B 2 fk^ 4 - 3p 



M <* )= tUJ • (30) 

Integrating the steady state magnetic energy spectrum from 
ki to the viscous wave number k v then gives 



M.i 



pk u I ki 

M{k)dk = B 2 I cr^da. 



The microphysical viscosity v must equal 

1 = (k^ v ' 

T(k v )k 2 v ki \ki 



v T (k v ) 



(31) 



(32) 



3.2 Invalidating the 3-D Zeldovich relations 
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where we have used dl9l and J2 II . This implies 



where 



Rm,i = Pr M-r — = Ptm t- 



p-3 



or 



Ra 



i.i 



Pr 



M 



l/(3-p) 



(33) 



(34) 



where Rm,1 = ^jrx~ is me magnetic Reynolds number asso- 
ciated with scale 

Let us now focus on the case Ptm = 1 to make con- 
tact with previous work. In this case, Em — Em,u = 
ll" M(k)dk. Then using ( 1341 in &3H gives, for p — 5/3 
(Kolmogorov spectrum), 

^B 2 lnR M A35) 



{b 2 ) = 2 / M{k)dk = 2B ln(k d jk x 

Jki 

For p / 5/3we instead find 



(l> 2 



2B 2 
5 — 3p 



5-3p 

t 1 1 



OR 2 5-3p 

5Zj**#. (36) 



for kd ki. Eqs. J35i and d36i are the generalized Zeldovich 
relations of Z83. 

Note that for p = 1 the result Eq. d36> would then imply 

(b 2 ) ~ iJA/^B 2 , (37) 

as noted by Z83, a formula consistent with the 2-D result d15t 
for L ~ I, However, as we discussed in section 2, there is no 
true steady state in 2-D and the derivation of i3bt relied on 
the steady state balance between the first two terms on the 
right of 1241 . Moreover, as also discussed below Eq. (9), the 
2-D turbulent kinetic energy spectrum is more consistent with 
p = 3, a value steeper than the Kolmogorov value p = 5/3, 
not more shallow. 

Before discussing the problems with the 3-D derivation 
of Z83 and the implications, we summarize a quite different 
derivation presented in KR that also gives a different result. 

3.1.2. The time dependent KR approach 

KR (p. 105; see also Low 1972) also start with (|23} with a 
fixed B, but instead of replacing the nonlinear terms by a dif- 
fusion term as in (1261 . they employ the quasilinear approx- 
imation (= second order correlation or first order smooth- 
ing approximation (FOSA)) to drop the nonlinear terms al- 
together. Furthermore, they keep the A dissipation term and 
the time evolution term on the left hand side of i26\ — both of 
which were dropped in the approach of the previous section. 
They then Fourier transform J26I in both space and time to 
obtain 

^ijk^klm B m ikj V[ 



bi = 



(38) 



— iuj + Xk 2 ' 
where the tilde indicates Fourier transform. Eq. d38l facili 
tates a rigorous relation between magnetic and velocity two 
point correlators. After some algebra, for isotropic, incom 
pressible turbulence this elegantly leads to (KR) 



M) 



k 2 Qu(k,u;) 



A 2 fc 4 



and Qij (k, uj)S(k + k', 



dkdu, (40) 

w') = (^(k,^-(k'V)}. 
limit, this gives 



(41) 



—2 



(O = f(R M )B 



(39) 



Defining s = cu/X and taking the A 

f(R M ) ~ — / Q u (k,0)dkds - -j- = R M - 

For A — > oo, KR instead obtain / = R 2 M . 

The KR approach produces a different (and spectra inde- 
pendent) result for f(Rnj) from the Z83 approach of the pre- 
vious subsection because KR keeps the time derivative term 
in (I23t and the resistive term, and drops the nonlinear terms. 
In contrast, Z83 drops the time derivative term and the re- 
sistive term and replaces the triple correlations by a turbulent 
diffusivity. If KR were to include the nonlinear terms by a tur- 
bulent diffusivity, the effect would be to replace A in this sub- 
section by a turbulent diffusion coefficient. The strong Rm 
dependence in \A\l would disappear, and the result would be 
spectrally dependent. 

3.2. Invalidating the 3-D Zeldovich relations 

Despite their differences, both approaches to deriving 3-D 
Zeldovich relations in the previous two subsections are miss- 
ing the same key ingredient which ultimate invalidates them 
for realistic steady states. Both approaches do not include the 
growth of the small scale field from the third term in (1231 . 

Fundamental to KR approach, was ignoring any contribu- 
tion from the third term in d23i via the second order correla- 
tion approximation or FOSA. But this excludes all the non- 
linear terms that make a nonlinear turbulent spectrum. These 
terms are not ignorable and as soon as b grows at all beyond 
the fixed B, they become dominant. 

While the Z83 approach of section 3.1.1 does not ignore 
the third term in d23l . it treats this term as a diffusion term, 
and incorrectly assumes that a steady-state arises between the 
first two terms on the right of (124b . This is the central prob- 
lem with the derivation of section 3.1.1 and invalidates d36i 
for realistic steady MHD turbulence. As discussed earlier, for 
steady state MHD turbulence initiated with a weak field, the 
small scale field can build up even in the absence of a large 
scale field. The saturated magnetic energy thus has nothing 
to do with a mean field, and Eq. d36l > is ruled out for the 
steady state. This discrepancy between d36i and numerical 
simulations of incompressible MHD turbulence (e.g. Maron 
& Blackman 2002; Haugen et al. 2003,2004; Maron et al. 
2004) arises because the third term of \2Al is a actually a 
source of both growth and decay (e.g. Kulsrud & Anderson 
1992; Haugen & Brandenburg 2004). The third term of OS) 
and should therefore be replaced by the sum of a growth term 
and a decay term. That is, 

(b ■ V x (v x b)) (fe) = Yb 2 {k) - iy T (k)k 2 b 2 (k), (42) 

where T ~ v(k)k indicates a growth rate. For a given k 3> 
k\, the terms on the right can largely balance independent of 
any term with the mean field or magnetic diffusivity. Accord- 
ingly, the energies of v and b could equilibrate without refer- 
ence to the second term in d24b which contains B. If the only 
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initial seed field were a weak mean field, then for a very short 
time, the second term of d24b would be the dominant term. 
For a k = mean field in a box with periodic boundaries, no 
change in the mean field is possible. Thus the second term of 
(1241 would swiftly become subdominant for an initially weak 
mean field, in contrast to that implied by the crucial step Eq. 
d26l , which appears in Z83. 

If instead of being a weak seed, the initial B were im- 
posed to be a fraction of order unity of equipartition with the 
turbulent kinetic energy, then the first term on the right of 
(I24> can be always competitive with the remaining terms on 
the right of (1241 . In this case, the small scale dynamo action 
is modestly suppressed by the strong mean field (Figs. 2 and 
6 of Haugen & Brandenburg 2004; also Brandeburg 1993). 

There is another subtlety in the derivation of (1361 . Ex- 
panding out the term containing B in ( 1241 gives 

(b • V x (v x B)) = (b q d n v q )B n - (b n v q )d q B n . (43) 

If correlations involving arbitrary products of small scale ve- 
locity and magnetic fluctuations were isotropic (that is, if the 
velocity and magnetic fluctuations were jointly distributed 
isotropically), both of the terms on the right would vanish: 
the first because the correlation is a vector, and the second 
because isotropy would make that term proportional to 5 qn . 
The value of the terms on the right hand side of d43i for non- 
helical flows and a weak mean field is therefore less than a 
simple order of magnitude estimate of those terms would in- 
dicate. Accordingly, the second term of J24i will be smaller 
than d25l > would indicate. Symbolically, Eq. d25i should be 
replaced by 

|(b • V x (v x B)) (fe) | ~ \b(k)kv(k)\ a B, (44) 

where the subscript a indicates that only an anisotropic part 
of the correlation contributes. This issue does not arise for the 
second term on the right of (I24i (which was taken to balance 
the left side of J44I in ( I28» which is given by 

(b • V x (v x b)} = (b q b n d n v q ) - (b n v q d q b n ). (45) 

Although the second term on the right of ( 145 \ vanishes for 
incompressible turbulence with the joint isotropy assumption, 
as can be seen by pulling the divergence out of the average, 
the first term on the right of ( 145 i has a surviving isotropic 
contribution. Thus, the fact that d43i might be even smaller 
than estimated in i25\ . further highlights that the terms on 
the right of J42I could largely maintain the steady state of 
i24\ by themselves. 

4. Conclusion 

We have re-examined the meaning and derivation of the Zel- 
dovich relations (Z57, Z8 3, ZR) of M HD for 2-D and 3-D. 
These relations (e.g. Eqs. ( I35I36I37I39I >) have been purported 
to relate the fluctuating magnetic field energy to a large scale 
or mean magnetic field energy in steady state forced MHD 
turbulence with an imposed mean field. We have identified 
why these relations do not apply for realistic steady states 
and highlight some common misconceptions about them. 

First, we have shown that these relations focus on the re- 
sponse of the total magnetic energy to non-helical turbulence 



driven by kinetic energy in the presence of an initially weak 
large scale field. Under the circumstances of their standard 
derivations, the large scale field serves only as a seed field, 
and does not grow. The relations say nothing about how large 
a mean field can grow with respect to a fluctuating field in he- 
lical MHD turbulence and therefore place no constraints on 
mean field helical MHD dynamo theory. 

Second, we point out that the presence of a mean field 
as a seed field in 3-D is the only reason that the small scale 
field that results from MHD turbulence would depend on it. 
Were the initial field a small scale seed field with no net 
mean, the small scale field would still grows, as seen in non- 
helically driven MHD turbulence simulations (e.g. Maron & 
Blackman 2002; Haugen et al. 2003,2004; Maron et al. 2004) 
Schekochihin et al. 2004), and the saturation value could not 
possibly depend on the mean field. For 3-D, we have indeed 
shown that standard derivations of the Zeldovich relation do 
not include (KR) or replace key nonlinear terms in the mag- 
netic energy equation by a turbulent diffusion term (Z83), 
ignoring an equally important growth term. The Zeldovich 
relations that emerge therefore result from an inappropriate 
balance of non-dominant terms. Accordingly, as seen in 3-D 
simulations, the ultimate steady-state saturation of the small 
scale field is determined by equipartition with the turbulent 
kinetic energy, independent of the mean field or magnetic 
Reynolds number for large Rm non-helical systems. 

Finally, for 2-D, we have shown that when the magnetic 
field has no k — component (but can still have a large scale 
k < kf component) the net magnetic energy must decay for 
a forced 2-D MHD turbulence system. Then the Zeldovich 
relation (I15t or (1371 cannot represent a true steady state. We 
have shown that instead, this relation represents the maxi- 
mum that the fluctuating field could obtain if the initial seed 

— 2 

field were a large scale seed field (and if RmB were less 
than the turbulent kinetic energy — otherwise the latter deter- 
mines a quasi-steady saturated value) after which the mag- 
netic energy would decay. In this respect, the import of the 
Zeldovich relation of Z57, was mainly to illustrate quantita- 
tively, the now well known result that 2-D MHD turbulence 
does not sustain a steady-state dynamo without a fixed seed. 
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